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ABSTRACT 

The geometric calibration of the Planck satellite using the planetary transits is in- 
vestigated, together with the reconstruction of any offsets from the nominal layout of 
the focal plane. The methods presented here may be applied to a single focal plane 
transit of a planet, to find the values of the geometric-calibration parameters at the 
epoch of the transit or all the transits over the course of the mission. The pointing 
requirements are easily met, with the pointing reconstruction being dominated by the 
errors due to the star tracker. 
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1 INTRODUCTION 

Planck is a European Space Agency satellite designed to 
produce high-resolution temperature and polarisation maps 
of the CMB. It possesses detectors sensitive to a wide range 
of frequencies from 30 to 857 GHz, split between two in- 
struments the HFI and LFI, the high and low frequency 
instruments, respectively. 

Planck will be inserted into a Lissajous orbit around 
the second Lagrange point of the Earth-Sun system, spin- 
ning about its axis once per minute. The line of sight being 
almost perpendicular to the spin axis, hence the detectors 
will almost follow a great circle. The spin axis will nomi- 
nally be repositioned every hour, and the roughly 60 or so 
circles corresponding to a single spin axis positioning may 
be binned together to form a ring. The nominal spin axis 
passes through the centre of the solar panels and is directed 
away from the sun, thus maintaining the rest of the satellite 
in a cone of shadow produced by the solar panels. The scan- 
ning strategy is determined by the sequence of the nominal 
positions of the spin axis over the course of the mission. 

The geometric calibration is the process whereby all the 
line-of-sight positions of the detectors are recovered, at any 
time during the observations. The relationship between the 
pointing and ti me may depend on multiple parameters, dis- 
cussed in full in Ivan Leeuwen et alJ (l2002t) : here we discuss 
only those parameters which require calibrating with the sci- 
ence data. The remaining geometric-calibration parameters 
may be determined solely from the star tracker data. 

The recovery of the geometric-calibration 
parameters has been investigated pr evious ly in 
lHarrison. van Leeuwen fc Rowan-Robinsonl j2004l) and 

* E-mail: dlh@ast.cam.ac.uk (DLH) 



lHarrison fc van Leeuwen] <l2005l) . which both use detec- 
tions from point sources in the ring data to recover the 
geometric-calibration parameters during and post-mission, 
respectively. In this paper we investigate the recovery of 
the geometric-calibration parameters using the solar system 
objects, specifically the planets. Additionally we show that 
it is possible to recover the focal plane layout, using the 
data corresponding to planetary transits. 

Section |2] contains an overview of the pointing re- 
construction of the satellite, in terms of the geometric- 
calibration parameters. Methods to reconstruct these pa- 
rameters and the focal plane layout from the detections due 
to the planets are presented in Section |3] The simulations 
used to assess the performance of these methods are de- 
scribed in Section 0] The results are presented in Section |S] 
where it is shown that the geometric-calibration parameters 
and any offsets from the nominal focal plane layout may be 
recovered using the transits of the planets. 



2 POINTING RECONSTRUCTION 

The pointing reconstruction of the satellite is achieved 
by the use of a star tracker. However, the exact rela- 
tionship between the reference frame which contains the 
star tracker and the satellite reference frame which con- 
tains the focal plane is uncertain. This relationship must 
therefore be established using the science data in order 
to meet the pointing accuracy requirements as dis cussed 
in lHarrison. van Leeuwen fc Rowan-Robinsonl i2004) . This 
may be achieved by the calibration of any offsets in the 
nominal values of the geometric-calibration parameters dis- 
cussed below. A more detailed discussion of the refer- 
ence systems and the attitude analysis may be found 
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Figure 1. The figure shows the geometric-calibration parameters 
which allow a description of the position of the field of view, 
FOV, with respect to the spin axis position, SA. The boresight 
angle, O-frp is the angle between the FRP and the spin axis. The 
rotation of the focal plane around the FRP with respect to the 
nominal scan direction is given by the roll angle, p. The reference 
phase, Vre/) is the value of the initial phase at the point at which 
the FRP crosses the reference point, as defined by the intersection 
of the reference circle and reference meridian. Where the initial 
phase is measured from the first point observed on the reference 
circle, and the reference point is given by the intersection point of 
the reference circle and the great circle connecting the spin axis 
position a the NEP. 



in Ivan Leeuwen et alJ i2002T) and lHarrison fc van Leeuwenl 
||200!tI) . 

The geometric-calibration parameters are defined to al- 
low the reconstruction of the lines-of-sight of the detectors 
with respect to the spin axis position, given the layout of 
the focal plane with respect to a reference point. This refer- 
ence point may be defined as the centre of the focal plane, 
and is henceforth referred to as the fiducial reference point, 
FRP. The geometric-calibration parameters may then be de- 
fined by three angles, shown in Figure Q The boresight an- 
gle is defined by the angle between the spin axis and the 
FRP, and the roll angle by the rotation of the focal plane 
around the FRP. The reference phase defines the position 
along the reference circle, and is measured from a refer- 
ence point defined by the intersection of the reference cir- 
cle with the great circle connecting the spin axis position 
and the North Ecliptic Pole, NEP. A more detailed discus- 
si qn of the geometric-calibr a tion p arameters may be found 
in lHarrison fc van Leeuwenl fcOOfJ) . 

If the focal plane layout is uncertain then as well as 
the geometric-calibration parameters, offsets in the individ- 
ual detector positions need to be recovered from the science 
data. These offsets may be expressed in terms of offsets in 
the nominal scan and cross-scan directions. 



3 METHODS 

The methods presented h e re are similar to those discussed in 
lHarrison fc van Leeuwenl |)200<t|) . the main differences arise 
from the fact that the positions of the planets are known to 
high precision. As previously, the geometric-calibration pa- 
rameters are assumed to be constant over the time frame of a 
ring and possess a slow linear variation over the course of the 
mission. The instantaneous offset in a geometric-calibration 
parameter, 7, may be defined on each ring by: 



7(r») 

where, 



U{Yi) = 



= 70 + 71 • U (r») 



(i) 



where Fi is the current ring number, Y max is the final ring 
of the mission and the ring numbers start from zero. The 
systematic offset in the parameter, 70, is hence defined as 
the instantaneous offset of the parameter exactly half-way 
though the mission and the drift in the value of the pa- 
rameter, 71, is defined as the total drift in the value of the 
parameter over the course of the mission. 

First the offsets in the reference phase and roll angle 
are extracted, then the offsets in the boresight angle. The 
geometric-calibration parameters values are then iterated 
over until convergence is reached, only then is the focal plane 
layout is recovered. 

3.1 Reference Phase and Roll Angle 

Offsets in the reference phase and roll angle may be found 
from the differences between the observed and expected 
phases of the detections, where the phase difference for a 
detection is given by: 



Aip = ip a —ip e 



(2) 



where ?p is the observed phase of the detection and i/> e is 
the expected phase for the detection. The expected phase 
of a detection on the i th ring in the d th detector may be 
calculated using: 



sin(/3 p ) — cos(ai) cos( 
sin(a 4 ) sm((pi) 



+ x d {pi) 



(3) 



where oti is the angular separation between the position of 
the planet, (X p , /3 P ), and the spin axis position, (Asa, /3sa), 
at the occasion of the i th ring, is the angle between the 
spin axis position and the NEP and Xd is the scan position 
of the detector relative to the FRP. 



The phase difference for a detection on the i 
the d th detector, Aipid, may also be given by: 



ring m 



At/Jid = i>ref + Ipreh ' k (T - y d {Pi) (PO + Pi ' ti (T i)) (4) 

where yd is the cross-scan position of the detector relative 
to the FRP at the time of the observation. The systematic 
offsets and drifts in the reference phase and roll angles may 
hence be solved for using a non-linear least squares analysis. 



3.2 The Boresight Angle 

The ordinate, r/, of a detection may be defined by the angular 
separation between the scan circle described by the detector 
and the position of the planet, and may be expressed as: 
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r\ = oti 



ad 



(•») 



where cti is the angular separation between the planet and 
spin axis positions and ay is the angle between the line-of- 
sight of the detector and the spin axis position and is given 
by: 



Old = Q.FRP (Ti) + Vd(pi) 



(6) 



where oifrp (fj) is the value of the boresight angle on the 



i th ring. 



By accumulating detections from multiple rings, and 
hence fitting the cross-scan transit of the planet through 
the detector, the ordinate corresponding to the peak of the 
cross-scan transit may be found. This corresponds to the off- 
set in the opening angle to this detector at the epoch of this 
planetary transit. The offset in the boresight is found for 
each planetary transit of the focal plane from the weighted 
mean of those off sets found for the detector s with the small- 
est beams, as in lHarrison fe van Leeuw cn (2005). In prac- 
tice this means the detectors which belong to the top four 
frequency channels. The systematic offset and drift in the 
boresight angle may be found by a straight line fit to the 
value found for each instantaneous offset from each focal 
plane transit. 



OlFRp(tp) = O>FRP + ClFRP 1 • tp 



(iv) 



(7) 



where T p is the mean ring number of the detections corre- 
sponding to those used in the evaluation of the boresight 
angle at the epoch of the planetary transit. 



3.3 Focal plane layout 

After the geometric-calibration parameters have been suc- 
cessfully recovered, the focal plane layout with respect to the 
FRP may be extracted. An uncertainty in the focal length 
of the optics will produce offsets in the focal plane layout 
which scale as a function of the detector position relative to 
the FRP. 

The offset in the position of a detector with respect 
to the FRP may be expressed in terms of an offset in the 
nominal scan, 5xd,o, and in the nominal cross-scan, Syd,o, 
directions. For each focal plane transit of a planet the offset 
in the detector position in the current scan and cross-scan 
directions are found. These may then be converted to offsets 
in the nominal scan and cross-scan direction using the value 
found for the roll angle at each epoch. 



Sxd,o = 8x d , p cos(p) - Sy d , p sm(p) 
5yd,o = Syd, P cos(p) + 6x d , P sin(p) 



(8) 



The systematic offsets in the position of a detector may 
then be found from a weighted mean of the offsets found, 
in the nominal scan and cross-scan directions, for each focal 
plane transit. 



4 SIMULATIONS 

In order to assess the performance of the methods devel- 
oped here in the reconstruction of the geometric-calibration 
parameters, it is necessary to simulate the planetary tran- 
sits occurring during the Planck mission. This requires an 
ephemeris and the brightness temperatures of the planets 



in the Planck frequency bands as well as a scanning strat- 
egy and information on the beams of the Planck detectors. 
Throughout this paper we have made the simplifying as- 
sumption of Gaussian beams. These methods, however, may 
be easily extended to non-Gaussian beams, due to detec- 
tor time constants or otherwise. Indeed, the beam profiles 
used are intimately connected to the geometric-calibration 
parameters, changes in the definition of the beam profiles 
will affect the geometric-calibration parameters. Hence, the 
geometric-calibration parameters are defined by the beam 
profiles used in their extraction. 



4.1 The scanning strategy 

While it is anticipated that this method will be applicable 
to any scanning strategy in which the circles corresponding 
to a single spin axis position may be binned together as a 
ring, only two potential scanning strategies for Planck were 
investigated here, a sinusoidal and a precessional scanning 
strategy. The sinusoidal scanning strategy may be described 
by: 



A fc = A + k9 
f3 k = Asm(n a X k ) 

and the precessional scanning strategy by: 
v = (A + kff) 



(9) 



sin(/3 fc ) = 
cos(<^>) = 

Afe = 
where, 

A() = A Q + 7T 



cos A 
cos(/3 fc ) 

v + - 
v — • 



' 2 



~ < n v v < ^ 



otherwise 



(10) 



(11) 



and Xq is the position of the sun at the time the first ring, 
k is the ring number, 9 is the angular separation between 
subsequent spin axis positions, n s , p is the number of periods 
within 2-7T, and A is the amplitude. The values of these pa- 
rameters used here are 8=2.5' , n a =2, n p = 2.05, and A=10° . 
This value of 0, given a repointing once per hour, keeps the 
spin axis directed away from the sun. 



4.2 The Planets 

An ephemeris code was written in order to evaluate the 
apparent positions and angular diameters of the planets 
as viewed from L2, the orbit of Planck about L2 is not 
included as this will not produce significant differences. 
The code uses the orbital elements available from JPL 
( http://ssd.jpl.nasa.gov/elem_planets.html) . 

Plan etary brightness te mpera tu res fromj d^Pater et^alJ 

ll989h. iGriffin fc Ortonl ll 19931). iGriffin et. al l 
Hildebrand et. al] lll985f). iLoewen stein et. all 
Muhleman fc Bergd ]l99lh. lOrton et. alT 




Rather Ulich fe Add 119741) . lUl ch fc Conklinl Jl976ft.lUlicrj 



Jl974 j) . IUlich. Coedell fc Davisl (Il973l) and lWhitcomb et. al 
( 1979) at far-infrared and sub-mm freque ncies were fitted 
using a third order polynomial fit, as in IGriffin fc Ortonl 
(1993), except in the case of Mars in which a straight line 
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Table 1. The brightness temperatures at each Planck frequency, 
used in the simulations to generate the observed fluxes for the 
detections due to the planets. 



Freq 


Mars 


Jupiter 


Saturn 


Uranus 


Neptune 


(GHz) 


(K) 


(K) 


(K) 


(K) 


(K) 


30 


211.7 


172.4 


143.9 


139.7 


135.8 


41 


212.0 


172.8 


144.7 


135.8 


131.8 


70 


212.5 


173.2 


146.0 


129.0 


124.8 


100 


213.1 


173.4 


146.9 


121.8 


117.5 


143 


213.9 


173.3 


147.4 


112.8 


108.3 


217 


215.4 


172.0 


146.2 


100.3 


95.6 


353 


218.0 


167.0 


139.0 


85.5 


80.6 


545 


221.8 


157.7 


124.0 


75.9 


71.7 


857 


228.0 


149.6 


112.7 


62.6 


62.2 



fit is sufficient. The brightness temperatures thus found for 
each planet, in each of the Planck frequency channels are 
shown in Table Q These brightness temperatures may be 
used to calculate the expected flux of each detection due to 
each planet, given the angular diameter of the planet at the 
time of the observation as provided by the ephemeris data. 



4.3 Generating simulated data 

As in Har rison fc van LeeuwerJ i2005li we directly simulate 
the list of detections, as delivered by the Level 2 DPC from 
their analysis of the time ordered data, TOD. This list of 
detections includes the position in phase and the amplitude 
observed for the transit, together with their respective er- 
rors. Additionally, the list includes the number of the de- 
tector which made the observation and the ring number on 
which it occurred. 

The simulations assume that the spin axis is reposi- 
tioned every hour and that the nutation effects are small 
so that the individual scans may indeed be combined into 
rings. For every spin axis position, the instantaneous val- 
ues of the geometric-calibration parameters are assessed and 
used to reconstruct the lines-of-sight of the focal plane. If 
a planet is located nearby, the amplitude with which it is 
observed is evaluated and if above a threshold signal-to- 
noise ratio it may be included in the simulated list of de- 
tections. Errors on the amplitude of the detections are gen- 
er ated assuming the goal values o f the noise in the ring, as 
in lHarrison fc van LeeuwerJ (|2005), and white noise. 

The magnitude of the error in the phase of the detec- 
tion, aj, , is assessed using the em pirical relationship found in 
lHarrison fc van LeeuwerJ i2005f) . However, this relationship 
must be modified to include the error in the velocity-phase 
relation, from the star tracker, which dominates for the high- 
est signal-to-noise detections. The value of is now given 
by: 

2 / 17 \ 2 2 2 

a + = \SNRj ab + a » (12) 

where SNR is the signal-to-noise ratio of the detection, at 
is the beam sigma of the detector and a v is the error in the 
velocity-phase relation, recovered from the star tracker data. 



Table 2. The errors in the offsets in the geometric-calibration 
parameters recovered at each of these planetary transits of the 
focal plane. 



Planet: 


Mars 


Jupiter 


Saturn 


Uranus 


Neptune 


Date : 


03/08 


04/08 


05/08 


06/08 


05/08 




(") 


(") 


(") 


(") 


(") 




0.06 


0.06 


0.07 


0.16 


0.28 


ffo FSP 


0.17 


0.23 


0.23 


0.38 


0.62 




3.71 


1.98 


3.51 


12.6 


25.2 



4-3.1 Saturation of detectors 

Jupiter may saturate the detectors; the effects of the po- 
tential saturation of the transits are considered here. Any 
saturation effects will only occur in extremely high signal- 
to-noise transits. The position of the peak of these transits in 
the scan direction is dominated by the error in the velocity- 
phase relation so the exclusion of the saturated measure- 
ments around the peak of the transit does not affect the 
accuracy at which the position of the peak of the transit 
may be found. 

As the planets are in the ecliptic plane the relatively 
large angular separations between successive rings in the 
cross-scan direction has the effect that only a small propor- 
tion of the observations in a cross-scan transit are saturated. 
If the saturated measurements are discarded then the posi- 
tion of the peak of the cross-scan transit may be successfully 
recovered. This again does not have a significant effect on 
the accuracy of the position of the peak of the transit. 

These effects were investigated using single scan and 
cross-scan transits; as the effects are not significant for any 
reasonable saturation limits they have not been included in 
the full simulations. 



5 RESULTS 

The simulated list of detections may then be analysed us- 
ing the methods presented in Section [3] Unless otherwise 
stated 500 noise realisations are performed for each simula- 
tion. The magnitude of the errors used are as follows, the 
velocity-phase relation, ov, and errors in the spin axis posi- 
tion coordinates, o-/s SA and &\ SA , are all assumed to be 1". 
The nutation amplitude used in the simulations was l"and 
the minimum signal-to-noise ratio of a detection used in the 
analysis is 30. 

Given the large number of high signal-to-noise ratio de- 
tections as a result of a planetary transit of the focal plane, 
and the high accuracy to which the positions of the planet 
are known it is possible to solve for the instantaneous values 
of the geometric-calibration parameters at the epoch of the 
focal plane transit. Table shows the errors in the recov- 
ered values of the geometric-calibration parameters at the 
epoch of a focal plane transit of each of the planets visible 
by Planck. Table |2] also shows that if there are only system- 
atic offsets in the geometric-calibration parameters then a 
single transit of any planet will reach the required pointing 
accuracy. The transits of Uranus and Neptune do not illumi- 
nate all the focal plane, which is the reason behind the large 
increase in the error in the roll angle. Due to the relative 
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Figure 2. These plots relate to the recovery of the line-of-sight, LOS, of a detector in the nominal scan and cross-scan directions, with 
the top panel in each plot containing the results for a 217 GHz, HFI detector and the bottom panel a 30 GHz LFI detector. Top left: the 
difference between the input and recovered value of the offset in the LOS of the detector in the nominal scan direction is plotted against 
the input value. The grey line is the global mean of these differences and the dashed lines enclose the region with ltr of this global mean. 
Top right: the difference between the input and recovered value of the offset in the LOS, in the cross-scan direction, is plotted against 
the input value. The grey line is the global mean of these differences and the dashed lines enclose the region with lrr of this global mean. 
The recovery of the offset in the cross-scan LOS, may hence be seen to be biased by 0.02" in the case of the 217 GHz detector and 
1.45" in the case of the 30 GHz detector. Bottom left: shows the dispersion in the mean recovered value of the offset in the nominal scan 
direction, black crosses, and the mean calculated error of the error in the recovered offset in the nominal scan direction, grey crosses, 
against the actual value of the scan offset. The dashed black lines enclose the ltr region about the mean dispersion, and the grey line 
is the mean value of the calculated error. Bottom right: shows the dispersion in the mean recovered value of the offset in the nominal 
cross-scan direction, black crosses, and the mean calculated error of the error in the recovered offset in the nominal cross-scan direction, 
grey crosses, against the actual value of the cross-scan offset. The dashed black lines enclose the 1<t region about the mean dispersion, 
and the grey line is the mean value of the calculated error. The calculated value of the error in the cross-scan LOS offset is seen to be 
an overestimate of the actual error in the recovered values. 



motions of Planck and Mars there are large variations in the 
time between successive transits of Mars, the first and po- 
tentially last opportunity to observe Mars transit the focal 
plane is March 2008. After which the next time Mars will be 
in the right relative location, will be in October 2009, after 
the end of the mission. 

Table shows the errors in the recovered values of the 
systematic offsets and the drifts in the geometric-calibration 
parameters using all the planetary focal plane transits over 
the course of the mission. Where the scanning strategy is as- 
sumed to start in July 2008, meaning no transit of Mars is in- 



cluded. Using an earlier start date for the scanning strategy, 
in order to include the March 2008 transit of Mars, does not 
however significantly improve the values displayed. Table |3] 
also shows that the scanning strategy employed has virtu- 
ally no impact on the ability of these methods to recover the 
geometric-calibration parameters, with no significant differ- 
ences found between the errors in the recovered values in the 
cases of the sinusoidal and precessional scanning strategies. 

Table 01 shows the effect of the errors in the geometric- 
calibration parameters in Table [3] on the reconstruction on 
the line-of-sight of the Planck detectors. The maximum error 
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Table 3. Comparing the errors in the recovered values of the 
geometric-calibration parameters in the cases of the sinusoidal 
and processional scanning strategies. In both cases only detec- 
tions with signal-to-noise ratios of 30 or greater were used in the 
analysis. 



Scanning Strategy: 


sinusoidal 


precessional 




(") 


(") 




0.03 


0.03 


°V>re/l 


0.13 


0.13 




0.11 


0.10 


eT °>FRP 1 


0.49 


0.47 


°"P0 


1.20 


1.19 


0> x 


4.82 


5.27 



Table 4. The maximum error in the line-of-sight of an HFI and 
LFI detector, due to the errors in the recovery of the geometric- 
calibration parameters shown in Table 151 



Scanning Strategy: 


sinusoidal 


precessional 




(") 


(") 


HFI (545 GHz) 


0.29 


0.28 


LFI (44 GHz) 


0.36 


0.36 



in the line-of-sight will occur at the beginning and end of the 
mission due to the higher sensitivity to errors in the drift in 
the geometric-calibration parameters. The location of the 
detector in relation to the FRP will also affect the error in 
its line-of-sight due to the uncertainties in the roll angle. 
Table [I] shows the largest errors found in the line-of-sight of 
the HFI and LFI detectors, corresponding to the values in 
Table 01 

Figure [5] demonstrates the ability of the methods pre- 
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Figure 3. The filled grey circles are the offsets in the detector 
positions. These offsets are 0.01% of the detectors position wrt 
to the FRP. The black crosses are the mean recovered positions 
found for the detectors with the lcr errors in the mean recovered 
value. There remain some small biases in the recovery of the de- 
tectors, especially in the nominal cross-scan direction of the order 
of a few hundredths of an arcsecond. 



sented here to recover the offsets in the line-of-sight of a 
detector. Offsets in the line-of-sight position of a detector 
with respect to the FRP are varied to assess the perfor- 
mance of the method. In each case the geometric-calibration 
parameters are also recovered, and are consistent with the 
zero offsets input to the simulations. The top panel in each 
plot, in Figure [5] corresponds to a 217 GHz HFI detector 
whereas the bottom panel corresponds to a 30 GHz LFI de- 
tector. The top left plot in Figure shows the difference 
between the input and recovered value of a line-of-sight off- 
set in the nominal scan direction against the value input to 
the simulations, whereas the top right plot shows the differ- 
ence between the input and recovered value of a line-of-sight 
offset in the nominal cross-scan direction against the value 
input to the simulations. In both plots the grey line shows 
the global mean of the differences and the dashed lines en- 
close the la region around the global mean. It may be seen 
in the top right plot in Figure |5] that the recovered value of 
the offset in the nominal cross-scan direction for each of the 
detectors is biased. This bias was found to depend on the 
number of planetary focal plane transits used in the analy- 
sis and the beam size of the detector, with the largest biases 
seen in detectors with the largest beams. It is this bias in 
the recovered cross-scan position which necessitates the lim- 
itation on the evaluation of the boresight angle to the top 
four frequency channels, which have the smallest beams. In 
the top right plot, the 217 GHz detector, FWHM of 5', is 
seen to be biased at the level of 0.02" whereas the 30 GHz 
detector, FWHM of 34.44', is found to be biased at the level 
of 1.45". Given that the magnitude of the bias is constant 
given the beam size and the planetary transits observed, it 
may be evaluated and removed if necessary, though at these 
magnitudes this is not required. The bottom left plot in Fig- 
ure |5| shows the errors in the recovered values of the offset 
in the nominal scan direction against the input offsets. This 
figure shows both the mean calculated values of the errors, 
grey crosses, and the dispersion in the recovered value of the 
offset, black crosses. The dashed lines enclose the la region 
around the global mean of the dispersion, and the grey line 
shows the global mean of the mean calculated error. The 
bottom right plot in Figure^ shows the errors in the recov- 
ered values of the offset in the nominal cross-scan direction, 
again both the dispersion in the recovered values and the 
mean calculated errors are shown. Here, however, the mean 
calculated errors are seen to be overestimated with respect 
to the dispersion in the recovered values. In both these fig- 
ures there is seen to be no dependence on the input value of 
the offset and the error in the recovered value. 

Since the bias in the recovery of the nominal cross-scan 
position of a detector is constant given the planetary transits 
observed, the bias may be removed. This is shown in Figure^] 
were the input offsets, filled grey circles, in the detector po- 
sitions are determined as a function of their distance from 
the FRP. These offsets are successfully recovered and any 
remaining biases are of the level of a few hundredths of an 
arcsecond. The black crosses show the recovered position of 
the detectors and the la error in that position. 

The total error in the line-of-sight of a detector includes 
the error due to its recovered position relative to the FRP. 
Table |5] shows the largest errors found in the line-of-sight of 
the HFI and LFI detectors, once the errors in the positions 
of the detectors relative to the FRP are included. 
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Table 5. The maximum error in the line-of-sight of an HFI and 
LFI detector, due to the errors in the recovery of the geometric- 
calibration parameters shown in Table 131 including the errors in 
the positions of the detectors with respect to the FRP. 





sinusoidal 


precessional 




cross-scan scan 


cross-scan scan 




(") (") 


(") (") 


HFI (217 GHz) 


0.34 0.22 


0.32 0.21 


LFI (30 GHz) 


0.43 0.41 


0.44 0.44 



6 DISCUSSION 

The geometric-calibration parameters and focal plane layout 
may be recovered with exquiste accuracy from the planetary 
transits of the focal plane. However, the transits of the focal 
plane by the planets are far from uniformly distributed in 
time. Depending on the start date of the mission there may 
a period of ~ 5 months without a planetary transist and it is 
possible that this period may occur at the end of the mission. 
If the geometric-calibration parameters vary with time then 
even if this is only a linear variation, it may be necessary to 
confirm this variation using the bright extra galactic point 
sources, as in lHarrison fc van Leeuwenl (120051) . An indepen- 
dent measurement of the geometric-calibration parameters 
may also be highly desirable in the case of an alteration to 
the scanning strategy, in order to improve the mapping of 
the beams of the detectors, when observing the planets 
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